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Abstract: Relativistic outflows carrying large scale magnetic fields have large inductive potential and 
may accelerate protons to ultra high energies. We discuss a novel scheme of Ultra-High Energy Cosmic 
Ray (UHECR) acceleration due to drifts in magnetized, cylindrically collimated, sheared jets of powerful 
active galaxies (with jet luminosity > 10 46 erg s" 1 ). A positively charged particle carried by such a 
plasma is in an unstable equilibrium if B ■ V x v < 0, so that kinetic drift along the velocity shear would 
lead to fast, regular energy gain. The highest rigidity (ratio of energy to charge) particles are accelerated 
most efficiently implying the dominance of light nuclei for energies above the ankle in our model: from a 
mixed population of pre-accelerated particle the drift mechanism picks up and boosts protons preferably. 



Introduction 

There is a consensus that relativistic outflows, and 
AGN jets in particular, are accelerated to relativis- 
tic speeds and collimated by large scale magnetic 
fields threading accretion disk and central black 
hole. Energetically, magnetic fields may carry a 
large fraction of jet luminosity [1, 3]. At largest 
scales magnetic field is dominated by a toroidal 
component B^. In addition, toroidal magnetic 
field may provide jet collimation, so that asymp- 
totically AGN jets may also be fully collimated to 
a cylindrical shape [2]. Axial motion of toroidal 
magnetic field create radial (in cylindrical coordi- 
nates) electric field . One also expects that jets are 
sheared, so that the central spine of the jet is mov- 
ing with larger velocity than its periphery. In a 
sheared jet Lorentz transformation cannot get rid 
of electric field everywhere in space. 

Powerful astrophysical outflows carrying large 
scale magnetic fields poses large inductive poten- 
tials that may be able to accelerate UHECRs [4] . If 
the Poynting luminosity of a source is Lem, then 
the total inductive potential is 
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Though Poynting luminosity L em may exceed the 
observed total luminosity of a source (which is re- 
lated to the rate of dissipation), the estimate (1) 
excludes acceleration of UHECRs in low power 
AGNs (e.g. Cen A and M87), low power BL Lacs 
and starburst galaxies (e.g. M82 & NGC 253) and 
limits the possibilities to more distant high power 
AGNs like higher power FR I, FR II radiogalaxies, 
radio-loud quasars and GRBs. 

Another constraint that acceleration sites should 
satisfy is that radiative losses should not de- 
grade particle energy. We can derive very gen- 
eral constraints on possible location of cosmic 
ray acceleration just by balancing the most efficient 
acceleration, by E ~ B, and radiative losses, f we 
normalize total energy density to energy density of 
magnetic field u — (tig, £ > 1, we find 
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where r c = e 2 /mc 2 . This gives limits on magnetic 
field and the distance from the central source 
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Where we also allowed a possibility that plasma 
is expanding with bulk Lorentz factor T, so that 
magnetic field in the plasma rest-frame is B /V and 
typical scale is R/T, and assumed £ = 1; EeV = 
10 18 eV. 

Relations (3) show that higher energy cosmic 
rays are better accelerated at larger distances. AGN 
jets, which propagate to more than 100 kpc dis- 
tances present an interesting possibility. Note, that 
as long as the jet remains relativistic, the total in- 
ductive potential is approximately conserved, so 
one can "wait" a long time for a particle to get ac- 
celerated without worrying about radiative loses. 
Thus, UHECRs can be accelerated inside the jet at 
distances from a fraction of a parsec (Eq. (3)) to 
hundreds of kpc, as long as the jet remains rela- 
tivistic and sustains a large inductive potential. 

Since inductive electric fields are orthogonal to 
magnetic field , particles cannot move freely along 
them. Thus, it is not obvious how to achieve en- 
ergy gain. One possibility is kinetic drift which 
may result in regular motion across magnetic 
field and along electric field leading to regular en- 
ergy gain as compared to the stochastic, Fermi- 
type schemes. Since the drift velocity increases 
with particle energy, the rate of energy gain will 
also increases with particle energy. Thus, the high- 
est energy particles will be accelerated most ef- 
ficiently. This means that from a pre-accelerated 
population, the mechanism proposed here will pick 
up particles with highest energy and will boost 
them to even higher values. In addition, when 
a particle has crossed a considerable fraction of 
the total available potential, it gyro-radius becomes 
comparable with flow scale. In this case, drift ap- 
proximation brakes down. As a result, as long as 
the particle remains inside a jet, the acceleration 
rate reaches the theoretical maximum of an inverse 
gyro-frequency. Thus, the most efficient accelera- 
tion occurs right before the particle leaves the jet. 
One may say that in case of inductive accelera- 
tion, becoming unbound is beneficial to accelera- 
tion, contrary to the case of stochastic acceleration 
when for unbound particles acceleration ceases. 



Particle dynamics in sheared flow 

In a transversely sheared flow one sign of charges 
is located at a maximum of electric potential, as we 
describe in this section. Consider sheared flow car- 
rying magnetic field . At each point there is electric 
field E = —v x B/c, so that the electric potential 
is determined by 

A$ = - V ■ (v x B) (4) 

c 

In a local rest frame, where v = 0, this becomes 
vanishes at the position of a particle, 

A$ = - (B • V x v) (5) 

c 

Thus, depending on the sign of the quantity 
(B • V x v) (which is a scalar) charges of one 
sign are near potential minimum, while those with 
the opposite sign are near potential maximum, see 
Fig. 1. Since electric field is perpendicular both 
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Figure 1 : Electric field in a sheared flow in a frame 
where surface the y — is at rest. For B • V x v < 
the electric field is directed away from the y = 
surface. 

to velocity and magnetic field, locally, the elec- 
tric potential is a function of only one coordinate 
along this direction. For (B ■ V x v) < (we will 
call this case negative shear) ions are near potential 
maximum. 

Inductive electric fields are not easily accessible 
for acceleration since particle need to move across 
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magnetic field , a processes prohibited under ideal 
Magneto-Hydrodynamics (MHD) approximation. 
On the other hand, kinetic drifts may result in reg- 
ular motion across magnetic field and along elec- 
tric field leading to regular energy gain. Typi- 
cally, the direction of drift is along the normal to 
the magnetic field and to the direction of the force 
that induces a drift. In sheared cylindrical jet with 
toroidal magnetic field , the electric field is in ra- 
dial direction, so that in order to gain energy par- 
ticle should experience radial drift. It is then re- 
quired that there should be a force along the axis. 
Such force may arise if a jet is axially inhomo- 
geneous, e.g. due to propagation of compressible 
waves along the jet, resulting in gradient drift due 
to changing magnetic pressure. 

Let's assume that there is a long wavelength iner- 
tial Alfven wave propagating along the z-direction 
with a phase speed Va ■ For Va <C c, the magnetic 
perturbation SB in the wave is much larger than 
electric perturbations SE by a factor c/Va, so that 
the wave is nearly magneto-static (in other words 
u <C k z c). A test particle will experience a drift in 
the x-direction with magnitude 
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where we assumed strong perturbation SB ~ B 
and uJsfl = eBo/mc. As a particle drifts along 
electric field in the ^-direction its Lorentz factor 
evolves according to 
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Thus, energy gain or loss of a test particle depends 
on direction of magnetic field, sign of charge 
(through iv b) and direction of velocity vorticity 
(through sign of Ly): in other words it depends on 
sign of shear. On the other hand, it is independent 
of the direction of the drift. 

Since drift velocity increases with particle energy, 
the rate of energy gain also increases with en- 
ergy, see Eq. (7). This leads to one of the most 
unusual properties of the proposed acceleration 
mechanisms: highest energy (or highest rigidity) 
particles are accelerated most efficiently. In addi- 
tion, at the last stages of acceleration, when parti- 
cle Larmor radius becomes of the order of jet scale, 



particle motion in positive shear flow becomes un- 
stable even without gradient drift while accelera- 
tion rate does reach absolute theoretical maximum 
of inverses relativistic gyro -frequency. 

In addition, since acceleration rate is inversely pro- 
portional to charge, at a given energy small charge 
(higher rigidity) particles are accelerated most effi- 
ciently. Thus, from a population of pre-accelerated 
particles with mixed composition, drift mechanism 
will pick up particles with smallest charge: pro- 
tons. This explains why above the ankle protons 
start to dominate over heavy nuclei. 

We can also calculate evolution of the spectrum. If 
the initial injection spectrum is power-law, f oc 
1 AS, then 
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Thus, for p > 2 the spectrum flattens with time. 

The hardest spectrum that can be achieved has a 
power law index of 2. This limiting case corre- 
sponds to unlimited acceleration in a plane-parallel 
geometry, which is realistically applicable to ener- 
gies well below the total available potential. At 
highest energies the final spectrum will depend on 
the distribution of pre-accelerated particles with re- 
spect to the electric potential and, in case of contri- 
bution from many sources, on distribution of total 
potentials. 

Discussion/predictions 

The best astrophysical location for operation of the 
proposed mechanism is cylindrically collimated, 
high power AGN jets. The proposed mechanism 
cannot work in spherically (or conically) expand- 
ing outflows since in this case a particle experi- 
ences polarization drift, which is a first order in 
Larmor radius, due to the fact that in the flow 
frame magnetic field decreases with time. For a 
constant velocity flow this drift is always against 
the electric field (for a positively charged parti- 
cle) and lead to the decrease of energy on time 
scale R/(cT), where R is a distance from the 
central source and V is the Lorentz factor of the 
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flow. Thus, in spherically expanding flows adia- 
batic losses always dominate over regular energy 
gain due to drift motion: the proposed mechanism 
would fail then. On theoretical grounds, AGN jets 
(or at least their cores) may indeed be asymptoti- 
cally cylindrically collimated [2]. Observations of 
large scale jets, e.g. Pictor A, do show jets that 
seem to be cylindrically collimated on scales of 
tens of kiloparsec. 

In cylindrically collimated parts of the jet acceler- 
ation can happen from sub-parsec to hundreds of 
kiloparsec scales: as long as the motion of the jet 
is relativistic the total electric potential remains ap- 
proximately constant. Thus, UHECRs need not be 
accelerated close to the central black hole where 
radiative losses are important. After a jet has prop- 
agated parsecs from the central source radiative 
losses become negligible. 

Another constraint on the mechanism comes from 
the requirement that in order to produce radial ki- 
netic drift the jet magnetic field should be inho- 
mogeneous along the axis. Though shocks pro- 
vide possible inhomogeneity of magnetic field (5- 
function inhomogeneity on the shock front), we 
disfavor shock since particles are advected down- 
stream and cannot drift large distances along shock 
surface. In a gradual inhomogeneity a particle 
drifts orthogonally to the field gradients and thus 
generally will remain in the region of inhomo- 
geneous fields. Extragalactic jets are expected 
to have axial inhomogeneities, both due to non- 
stationary conditions at the source and due to prop- 
agation of compression and rarefaction waves gen- 
erated at the jet boundary via interaction with sur- 
rounding plasma. 

Our model has a number of clear predictions, some 
of which are related to astrophysical association of 
acceleration sites of UHECRs with AGN jets [5] 
and some are specific to the model: (i) one needs 
a relatively powerful AGN, with luminosity > 
10 46 erg/sec. This limits possible sources to high 
power sources like FR II radiogalaxies, radio loud 
quasars and high power BL Lacs (flat spectrum ra- 
dio quasars). Powerful AGNs are relatively rare 
and far apart, so that a steep GZK cut-off corre- 
sponding to large source separation should be seen, 
(ii) UHECRs come from sources with low spacial 
density. This may be reflected in the distribution 
of arrival directions, (iii) Extragalactic UHECRs 



should be dominated by protons, (iv) Depend- 
ing on "extra-galactic seeing conditions" arrival di- 
rections of UHECRs may point to their sources, 
though complicated magnetic field structure may 
erase this correlation. In addition, the fact that 
only flows with negative shear can accelerate pro- 
tons implies that only approximately half of such 
AGNs can be sources of UHECRs (this assumes 
that the AGN central engine - black hole or an ac- 
cretion disk - is dominated by large scale, dipolar- 
like magnetic field). 
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